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The cumulants of the logarithm of the conductance (Ing) in the localized regime in the one- 
dimensional Anderson model are calculated exactly in the second Born approximation for weak 
disorder. Only the first two cumulants turn out to ne non-zero since the third and fourth cumulants 
vanish identically and the higher cumulants are of higher order in the disorder. The variance and the 
mean of Ing vary linearly with length L while their ratio is proportional to the inverse localization 
length. The resulting exact log-normal distribution of conductance thus corresponds to a special 
form of two-parameter scaling. This contradicts the standard one-parameter scaling in the random 
phase approximation. 
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The appearance of the scalins theory of localization in d-dimensional disordered systems [l[ and more detailed 
developments of it in ID |^ |M ^ |^ and for quasi ID systems [EH has inaugurated a golden age for mesoscopic 
physics. This is demonstrated by the numerous review articles Isffg UlOl llll Il^ll3| and monographs deaHng 
with application of scaling ideas to transport in disordered conductors. 

The fundamental hypothesis in the scaling theory P] is that the scaling of the logarithm of a typical conductance 
g as a function of a characteristic size L of the system is described asymptotically for large i by a universal function, 
f3(\ng), of a single parameter (SPS), namely In g itself: 

The function (3 which may generally depend on dimensionality is independent of L and of any microscopic parameter 
in the system. We recall that in the studies of scahng in ID systems j^i^ AiSj] the parameter ln((l -I- p), with P = 
the resistance, was identified as the convenient scaling variable both in the low resistance (p << 1) (quasi- metallic) 
regime and in the large resistance or insulating regime {p » 1). This variable reduces to — Ing in Eq. JQ) for p >> 1 
(localized regime) and, thanks to the Landauer formula, p = (with r^, and tj^ the reflection and transmission 
coefficients od the system, respectively), it coincides with — Int/^ 2± —Ing. 

In Refs. 0, IMS 13 on ID systems it was argued that the scaling theory of Abrahams et al. had to be interpreted 
in terms of the scaling of the distribution Pg{g) of the random conductance of the system. SPS then means that 
Pg{g) is fully determined by a single parameter such as e.g. the mean logarithm, (—Ing), which is itself defined by a 
scaling equation of the form with In g replaced by (Ino) . 

The validity of SPS in the theory of Abrahams et al. jl| has generated debates and controversies in the past, which 
have not been fully settled and which have recently been revived 0, 0, 0|. Indeed the justifications of the SPS 
hypothesis in the analyses |M S 13 rests on the use of a random phase approximation (RPA) which assumes that 
the phases of the amplitude reflection- and transmission coefficients Rl and Tl (with = |i?L|^,iL — |7lP) are 
uniformly distributed over (0, 27r)in the localization domain, which corresponds to length scales L much larger than 
the localization length ^. Although strong evidence, both numerical ^^'^ analytical 113,01, exists for uniform 
phase distributions for L >> ^ in the ID Anderson model [i^l, the controversy about the validity of SPS does not 
seem to be satisfactorily resolved 0,0, E,!!^- However, Cohen et al. have concluded that if RPA was not 
used one would obtain at most two-parameter scaling. Doubts about the RPA results have led Deych et al. 0| 
to reconsider the scaling problem in the framework of an exact solution for the Lloyd model, which differs from the 
Anderson model by the use of a Cauchy distribution of site potentials. While the study of the Lloyd model does 
not allow to draw conclusions for scaling in the Anderson model (because the Cauchy distribution possesses infinite 
moments), the results of Ref. 0| do show that the phase randomness model is invalid for the Lloyd model. On the 
other hand, Schomerus and Titov have studied the deviations of the conductance distribution from the SPS form 
obtained in the RPA at the band centre {E = 0). They also discussed using large deviation statistics, the 

higher cumulants which, beyond the mean and the variance, affect the exact distribution of conductance. 

The purpose of the present paper is to calculate exactly the distribution of conductance in the localization regime 
(L » ^) in the current ID Anderson model for weak disorder. We feel that, in particular, our results definitely 
clarify the status of the SPS hypothesis in the Anderson model. 
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The Schrodinger equation for a chain of N disordered sites 1 < m < of spacing a = 1 {L = N) is 

fn+l + ^n-l + eny^n ^ E (fn , (2) 

where the site energies e„, in units of a constant hopping rate, are independent gaussian variables of zero mean and 
correlation 

{em£n) = r]'^Sm,n ■ (3) 

The disordered chain is connected as usual at both ends to semi-infinite non-disordered chains — 0) with sites 
m > N and m < 1, respectively. 

The distribution of the transmission coefficient tjy (conductance) for an electron incident from the right with 
wavenumber —k (or energy E — 2cosfc) will be obtained by using the general recursion relations of Rcf. 21] which 
connect the transmission- (reflection-) amplitudes of a chain of n sites with the corresponding amplitudes for a chain 
with one less disordered site of length n — 1. By iterating the transmission amplitude recursion relation (11. a) of [ij 
we readily find 

JV 

lnrAr = zfciV-^ln[l-i//„(l-|-e2''=i?„_i)] , (4) 

n=l 

where the reflection amplitudes are given by Eq (11) of 



P _ e^'fefl„_i+z;/„(l + e^'fei?n_i 



with 



(6) 



2 sin k 

The Eq. Q is our starting point for studying the cumulants of In In = In Tn + In and its probability distribution. 

For weak disorder (small £„) we restrict the analysis to low order in the correlation parameter 77^. Since the variance 
of In i^v is found to be proportional to ry** to lowest order an exact treatment for weak disorder must involve expanding 
Q to 4th order in the site energies. In parallel our analysis requires the reflection amplitudes in Q to second order 
in the disorder, that is in the second Born approximation of backscattering. Since Rm is linear in the site energies to 
lowest we expand it in terms of linear and quadratic contributions, Rm — Rm'^ + Rm\ which are defined by recursion 
relations obtained by identifying contributions of the same order in an expansion of both sides of lO): 

Rl'^=e'"'Ri'l,+ii.n , (7) 



R^^^e^''R':\ + 2^.„e'^'^Ri'l,-.l . (8) 

Solutions of 1)718(1 will be used for extracting explicit forms of successive order terms up to fourth order in the expansion 
of In tpf obtained from Q , 

N 4 p 

InT^ = V V(-l)f+i^ , u„ = -zi.„(l -I- e'^'^R,,^,) . (9) 
^ — ' ^ — ' V 

n=lp=l ^ 

At this point we have to determine the relevant asymptotic solutions of (|7I8|I for describing transport in the strong 
localization regime, N » ^. In this domain the transport is dominated by the amplitudes i?„ at large scales n. 
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which are described by some invariant distribution y, which is independent of the initial values at sites where the 
iteration of lO was started. On the other hand, the probability densities associated with i?„ and Rn-i, respectively, 
are related by a linear integral equation defined from the recursion relation 0. Since for large n these densities are 
invariant it follows that in the recursion relation (jSJ we may replace by i?„ in this limit. Then we obtain from 

(TTIHl 



= (10) 



and 



R^?'-T^[j^--n] ,n»^ . (11) 



The existence of an invariant density for large n is a general property which is equivalent to the random Oseledec 
theorem (a review on these aspects and of their application in weak disorder studies is found in j25j|, see also 26]). 

Using 1)911 l|l we express the powers of \ntN , (IntAr)" , n = 1,2,3,4 to 4th order in the site energies (powers 
higher than 4 have only contributions of higher order in the site energies) in order to find the moments, using (jS)). In 
performing the averages we make use of the factorization property 



valid for arbitrary indices m,n,p,q different or not j27l |. A useful simplification is also to note that averages of the 
form {vnRn~i) vanish because depends only on energies of sites m < n — \. By straightforward but somewhat 
cumbersome calculations we obtain the following final results for the moments: 



(Inijv) = -27iV , 7 = ^^-3— fl ^1 , (13) 

^ ' Ssin^/fc 64sin4fcV 2sin2fc/ ^ ' 

{{\ntNf) = -^^{N^ + m), N»^ , (14) 
16 sm k 

((Intw)") =0(7;'^) , n>3 , (15) 

which are exact within the second Born approximation. The vanishing of the n = Z and n ^ A moments (Eq. (|15|) ) 
through order rj'^ results from the fact that the fourth order terms in both (Iu^at)^ and (Ini^r)^ cancel identically. 
Finally from 1)12(1 and H13|l we get for the second cumulant of In i^r : 

2n^N 

varlntAT^ \ = 127^iV + 0(77^) , (16) 
16 sm k 

and higher cumulants defined e.g. in vanish to fourth order, of course. 

In discussing the above results we first note that 7 = ■C""^ in l(T^ is nothing but the perturbation expression for 
the inverse localization length, in which the leading term coincides with the well-known expression in the (first) Born 
approximation j^^. Our results for the cumulants of Int/, ~ In 17 in the localized regime {N >> show that the 
distribution of the conductance has an exact log-normal form in the second Born approximation, namely 

1 (In g + 2-yJV)^ 

^^9)dg = ^==e ding . (17) 

V247r72iV 

The important new feature of the exact result for the variance in H16|l is that while being linear in iV, in agreement 
with the central limit theorem, the ratio of it to the mean ()13|l is proportional to the inverse localization length. In 
this case the criterion for SPS which is obeyed in RPA, is clearly violated. This shows that the mean and the 
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variance of In g in the distribution (|16(l act as distinct parameters even though they obey the same scahng equation 
of the form -j^^ = / j / = (l^iff) oi' varlngi. In this sense our exact analysis demonstrates two-parameter scahng in 
the locahzed regime of the Anderson model, in contrast to the RPA which leads to a strictly SPS H,IEQ>I1|' 

Our analytical results for the cumulants of Ing differ qualitatively from the results for the ID Anderson model 
obtained by Schomerus and Titov (ST) ^^Ij '^ho did not use RPA. ST find that at the band center {E = 0) the first 
three cumulants differ from zero and are proportional to (In g) , which indicates that the distribution of In g obeys a 
SPS. Furthermore, for large energies ST recover the SPS results of RPA ITj. We believe that the differences between 
our results and those of jl7l | originate in the fact that we are using an exact analysis of the discrete Anderson model 
whereas the treatment of ST deals with a continuum limit of the tight-binding equations. In particular, the disorder 
enters in our analysis via the reduced site energies ^ depending on the Bloch wavenumber k at energy E = 2cosfc. 
On the other hand, the length scale dependence of successive moments in (I13I15II enters, quite distinctly of the disorder 
parameters £«, via multiple summations over sites of the form (@J. In contrast, in the Fokker- Planck type continuum 
treatment of ST the disorder enters exclusively via a rescaled length and a reduced continuum energy, which causes 
the cumulants at = to depend on the rescaled length only (via (Ing)). 

Finally, we note that the observed shortcoming of the RPA conductance (resistance) distribution in the localized 
regime exists also in the well-known DMPK equation for the distribution of transmission eigenvalues in quasi ID 
systems This is because the DMPK equation reduces in the ID limit to results obtained in the RPA. 
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